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ABSTRACT 

We investigate HST/ ACS and WFPC2 images at the positions of five candidate microlensing 
events from a large survey of variability in M31 (MEGA). Three closely match unresolved sources, 
and two produce only flux upper limits. All are confined to regions of the color-magnitude diagram 
where stellar variability is unlikely to be easily confused with microlensing. Red variable stars 
cannot explain these events (although background supernova are possible for two). If these lenses 
arise in M31's halo, they are due to masses 0.15 < m/MQ < 0.49 (95% certainty, for a (5-function 
mass distribution), brown dwarfs for disk lenses, and stellar masses for bulge lenses. 

Subject headings: gravitational lensing — galaxies; individual (M31) — galaxies: halos — dark matter 



1. Introduction 

Galaxian dark matter has been recognized for 
over 70 years (Zwicky 1933), and tied in part to 
the halo for over 30 (Rubin & Ford 1970). Halo 
dark matter's nature is still a mystery. Gravi- 
tational microlensing can reveal individual dark 
matter objects of roughly stellar mass (Paczyhski 
1986). To test this, MACHO observed the Mag- 
ellanic Clouds for 5.7 years, (Alcock et al. 2000) 
and EROS (Afonso et al. 2003) did so for 5. The 
former report microlensing events more common 
than the known, purely stellar expectation, with 
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lensing fraction / 20% of the dark matter halo 
mass (8-50%, with 95% confidence) of -0.4 M© 
masses. EROS found / consistent with zero (but 
marginally consistent with / w 20%). 

M31 microlensing could potentially settle this 
quandary definitively (Crotts 1992). Since we 
can explore microlensing across the face of M31, 
we can use this distribution to distinguish where 
in the galaxy the lenses arise. Several surveys 
of M31 microlensing (Riffeser et al. 2003, Joshi 
et al. 2004, Calchi-Novati et al. 2005, including 
MEGA: de Jong et al. 2004 and its predecessor 
VATT-Columbia: Uglesich et al. 2004), together 
report ~20 probable microlensing events, and have 
a tendency to confirm the MACHO result. 

With its crowded target stars, M31 microlens- 
ing relies on image subtraction to reveal event 
lightcurves, which removes the baseline flux. Us- 
ing HST to recover the source flux (e.g. Ansari et 
al. 1999, Auriere et al. 2001), one can compute 
event amplification, hence Einstein parameters, 
constraining physical parameters e.g., lens mass. 
MEGA and VATT-Columbia also use source star 
color to distinguish microlensing from variable 
stars, since very red variables (miras and semireg- 
ulars) produce outbursts that, with their base- 
lines subtracted, mimic point-source, point-lens 
( "Paczynski" ) light curves (Uglesich et al. 2004). 
Residual flux from these events, however, is red- 
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dcr than almost all potential microlensing source 
stars. MEGA will soon publish its microlensing 
sample, and now is an excellent opportunity to 
check these event selection criteria. 

2. Observations and analysis 

To study candidate events we appeal to superior 
HST angular resolution: 160 ACS and WFPC2 
images, taken in F555W and F814W filters in 16 
orbits, cover 0.17 deg^, ~30% of the MEGA field. 
Here we study the largest sample of candidate 
microlensing events, the INT/WFC subsample of 
MEGA (de Jong et al. 2005), in order to under- 
stand and improve ground-based selection criteria. 

The analysis used is detailed by Cscrcsnjcs 
et al. (in preparation). We carefully align the 
HST and ground-based images by matching cat- 
alogs of ground-based versus Gaussian-convolved 
HST sources for each filter combination {HST, 
F555W/F814W versus INT, r'/i' and KPNO 4m, 
R/I), providing up to 8 different position esti- 
mates. For a given ground-based position, the 
two independent HST positions (via F555W and 
F814W) always agree to < 0".03 (typically 0".01), 
so positional accuracy depends mainly on ground- 
based data. The adopted position is a weighted av- 
erage of individual estimates (Fig. 1). The spread 
of different estimates for each event are 0".02 to 
0".08. Of five microlensing candidates analyzed, 
we identify three sources and find flux upper lim- 
its for two. 

HST data were photometered with DAOPHOT 

(Stetson 1987), as prescribed in Sirianni ct al. 
(2005). The locations of the candidate microlens- 
ing sources on a color-magnitude diagram are 
shown in Fig. 2. For each candidate event, we 
normalized the differential light curves to R-fluxes, 
using color-magnitude diagrams and HST base- 
line fluxes (for the two undetected events, using 
the baseline flux upper limit), then performed a 
Paczynski fit in {uo,to,tE) to the combined light 
curve (Fig. 1). For the two undetected events, the 
resulting tj^ corresponds to a lower limit. 

Only the Einstein time-scale (f^) constrains 
lens characteristics, particularly its mass to. For a 
given time-scale ts, the lens-mass probability dis- 
tribution is P{M,tE) = (dr/dt£;)/r, where T is 
the event rate (Gricst 1991). We consider alterna- 
tively a lens located in the halo of M31, in the disk 



or in the bulge. For halo lensing, we consider the 
simple case where the lens is part of a spherical 
isothermal halo composed of single mass objects, 
with a density distribution defined as 

where R is the radial distance to the center of 
M31, and Tc is a core radius of 5 kpc. The one- 
dimensional velocity dispersion of the lenses is set 
to 170 km s"'^, consistent with a rotation curve 
of 240 km s~^. The disk is modelled by a double 
exponential with a radial scale length of 6 kpc and 
a vertical scale length of 400 pc. The bulge corre- 
sponds to the "small" bulge model of Kent (1989), 
with a velocity dispersion of 150 km s~^. More de- 
tails about the model can be found in Baltz, Gyuk, 
& Crotts (2003). 

3. Individual events 

ML-8: this event's position lands within the 
FWHM of a red clump star's image, with R — I 
in excellent agreement with the peak flux's color 
in differential light curves (0.60 ± 0.16 vs. 0.59 
mag). With the baseline set to this star's flux, 
a Paczynski flt yields ampliflcation A = 8.49, an 
Einstein time-scale = 60.6 ± 4.2 days. The cor- 
responding lens masses are m = O.ntoilMQ for 
a halo lens, m = O.OStg^os-^© ^ *iisk lens, and 
TO = 2.85^2:3? Afg for a bulge lens. 

However, this event lands ^ 0".9 from the cen- 
ter of a background galaxy (subtending ~ 1".5 x 
0".3). Its color, flux and decline rate are consis- 
tent with a Type la supernova at z ?a 0.5, with 
< 1 mag extinction (see Johnson & Crotts 2005). 
One must balance the number of SNe (~ 100 y~^ 
deg~^ e.g.. Woods & Loeb 1998) landing within 
the FWHM disk of a source star of consistent color 
we would detect (1.2 arcsec"^) versus the number 
of microlensing events (evidently ^5) landing so 
close to an i? < 23 galaxy (100 arcmin^^ — Huang 
et al. 2001). The expected number of both kind 
of events are of the order of a few tenths, with 
perhaps microlensing being slightly more likely. 

ML- 10: this event lands within the FWHM disk 
of a giant branch star of color R — I = 1.05, in per- 
fect agreement with the microlensing data. This 
source has A = 4.00 and = 64.7 ± 1.9 days, cor- 
responding to a halo lens mass to = 0.33^0 23-^©' 
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Fig. 1. — Left to right: 1) (1".5)^ HST image (f814w band) around microlensing candidates. The circles 
correspond to la spread of the individual position estimates. The grid represents the INT pixel sampling; 
for ML-16, the cross point corresponds to an independent estimate as described in the text. 2) full combined 
light curve (Filled squares: KP-R, asterisks: KP-I, open circles: INT-r', open triangles: INT-i'); 3) Zoom on 
the event peak; 4) Lens-mass probability distribution for a lens in an isothermal halo (full line), in the disk 
(dashed line), and in the bulge (dotted line). 
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Table 1 

Event source photometry and Microlensing parameters 



id 


ML-8 


ML-10 


ML-11 (S4) 


ML-16 (Nl) 


ML-18 


R.A. (J2000) 


00:43:24.53 


00:43:54.87 


00:42:29.90 


00:42:51.22 


00:43:17.27 


Declination (J2000) 


41:37:50.4 


41:10:33.3 


40:53:45.6 


41:23:55.3 


41:02:13.7 




24.94 ± 0.14 


23.36 ± 0.09 


24.86 ± 0.30 


> 23.86 


> 25.09 


^hst 


24.34 ± 0.08 


22.31 ±0.07 


24.71 ±0.26 


>23.32 


>24.59 


{R-i)\hst 


0.60 


1.05 


0.15 






(R-i)k= 


0.59 


1.05 


0.21 




0.51 


-^max 


8.49 


4.00 


41.93 


>16.01 


>11.42 


t„/day 


60.6 ± 4.2 


64.7 ± 1.9 


26.1 ± 1.1 


>6.9 


>86.6 




0.89 


1.26 


1.01 


1.29 


1.04 


"Ihalo/M©"- 


^■•''--0.21 


U.OO_Q 23 




> 0.00 


> 0.62 


"»di8k/M0* 




0.05;°:™ 


01+" °^ 


> 0.03 


> 0.09 


"»bulge/M0* 


2.85;-°? 


2.oo;i:f, 


U.JO_Q 28 


> 0.05 


> 3.73 


comment 


red clump or SN 


giant branch 


very blue 


undetected(?) 


in cluster or galaxy 



^most likely lens mass (with la confidence intervals) 




ML- IS. 



ML-IB 



Fig. 2. - Color-magnitude plot of 5 candidate 
sources, with upper flux limits as thick lines. The 
dashed lines enclose the area where LPVs and 
semi-regulars are expected (Brown et al. 2004). 



a disk lens mass m = O.O5^q q4M0, or a bulge 
lens mass m = 2.00^\'ffMQ. It lands suspiciously 
close to a region of the CMD common to variables. 
Still, the achromaticity of the variation, the well- 
fit and well-sampled peak {x^ /N = 1.26), and the 
stability of the baseline over 7 seasons strongly in- 
dicate a real microlensing event. 

ML-11: this event lands on a faint blue star 
(R — I = 0.15 ± 0.40) severely blended with a red 
clump star. The light curves fit yields a similar 
R - 1 = 0.21. Its baseline flux implies A = 41.93 
and = 26.1 ±1.1 days. This event, from Paulin- 
Hcnriksson et al. (2002), is near M32, suggesting 
that the lens resides there. If not, the most likely 
lens mass is m 

1+0.02 



0-05j:^ o^M0 in the halo, m 

+0.31 : 



O.OlI'o:oi^0 in tlie disk, and m = 0.361^-^8-^0 in 
the bulge. 

ML-16: this event lands in a WFPC2 field and 
was also seen by POINT-AGAPE (Auriere et al. 
2001). They publish a color for the event peak 
based on INT g' and r' (no i' data are available), 
corresponding to V-I ^2.1. We find no detected 
source at this position; the nearest detected star 
landing ~ 0".l away (one WFPC2 pixel) with 
V — I ~ 1.1. Using the flux of this star as an 
upper limit, we find A > 16.01, > 6.9 days, 
and m > O.OOSM© being poorly constrained for a 
halo lens. 

Auriere et al. (2001) seem to have isolated a 
different source star. Their celestial coordinates 
disagree with ours by 3", but cannot be checked 
from published data since no image of the source 
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field is provided. In order to check our astrometry, 
we repeated the procedure we applied to r' data to 
g' data retrieved from the INT archive. We also 
repeated the same procedure (for both r' and g' 
data) using the WFPC2 images taken by Auriere 
et al from the archive. Finally, one of us (A.C.) 
made an independent check by choosing 8 bright 
unsaturated, isolated stars as coordinate inputs to 
IRAF geomap to construct the coordinate trans- 
form. These various estimates agree to better than 
0.5 pixel in the INT data, whereas the nearest star 
(two of them actually) with consistent colors and 
magnitudes to that claimed by Auriere et al. are 
0".6 (or almost two INT pixels) away. Our posi- 
tion is marginally consistent with the faint source 
cited above and in Table 1, but inconsistent with 
that of Auriere et al. 

ML-18: this event lands in a bright region, per- 
haps a cluster or background galaxy. We isolate no 
source here, so provide only an upper limit base- 
line flux, estimated by taking the brightest pixel 
// 

within .05 and considering that it contains at 

most 15% of the source flux, as constrained by 
the point-spread function for ACS. With this flux 
limit, A > 11.42, > 86.6 days, and m > 0.62 M© 
for a halo lens, m > 0.09 Mq for a disk lens, and 
m>3.73 Mq for a bulge lens. 

4. Conclusions 

Of five events in our fields, we find three likely 
matches, and baseline flux upper limits on the 
other two. Colors of the three identified sources 
agree with those obtained from their differential 
light curves alone. The two upper limits displace 
these events from the asymtotic giant branch, 
where confusing mira and scmircgular variables 
can occur. No candidate is a bright red variable. 
One might interpret ML-8 as a supernova, but a 
microlcnsing event is just as probable. We also 
cannot rule out a supernova as the source for ML- 
18, which might also coincide with a background 
galaxy. In a future paper, the complete MEGA 
data set will fill out ML-18's light curve; unfortu- 
nately, we have no additional data on ML-8. 

Taking the product of the individual mass prob- 
ability distributions obtained for each event, these 
lenses in a halo model (Baltz, Gyuk & Crotts 
2003) of a single component mass are constrained 
to 0.15 < m/MQ < 0.49 at the 95% level. M31 



microlcnsing rates may be consistent with pure 
self-lensing (de Jong et al. 2005), so we consider 
bulge lenses (0.64 < m/M© < 2.02), or disk lenses 
which correspond to probably unrealistic brown 
dwarf masses (0.02 < m/M© < 0.06). 

We acknowledge support from STScI (GO 
10273) and NSF (grants 0406970 and 0070882). 
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